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ABSTRACT 

Context. Extremely reddened asymptotic giant branch stars (AGB) lose mass at high rates of > 1(T 5 M 0 yr _1 . This is the very last 
stage of AGB evolution, in which stars in the mass range ~2.0-4.0 M 0 (for solar metallicity) should have been converted to C stars 
already. The extremely reddened AGB stars in the Galactic bulge are however predominantly O-rich, implying that they might be 
either low-mass stars or stars at the upper end of the AGB mass range. 

Aims. Our goal is to determine the mass range of the most reddened AGB stars in the Galactic bulge. 

Methods. Using Virtual Observatory tools, we constructed spectral energy distributions of a sample of 37 evolved stars in the Galactic 
bulge with extremely red IRAS colours. We fitted DUSTY models to the observational data to infer the bolometric fluxes. Applying 
individual corrections for interstellar extinction and adopting a common distance, we determined luminosities and mass-loss rates, 
and inferred the progenitor mass range from comparisons with AGB evolutionary models. 

Results. The observed spectral energy distributions are consistent with a classification as reddened AGB stars, except for two stars, 
which are proto-planetary nebula candidates. For the AGB stars, we found luminosities in the range ~ 3000 - 30,000 L Q and mass-loss 
rates ~ 1(T 5 - 3 x 1(U 4 M 0 yr -1 . The corresponding mass range is ~ 1.1 -6.0 M 0 assuming solar metallicity. 

Conclusions. Contrary to the predictions of the evolutionary models, the luminosity distribution is continuous, with many O-rich 
AGB stars in the mass range in which they should have been converted into C stars already. We suspect that bulge AGB stars have 
higher than solar metallicity and therefore may avoid the conversion to C-rich. The presence of low-mass stars in the sample shows 
that their termination of the AGB evolution also occurs during a final phase of very high mass-loss rate, leading to optically thick 
circumstellar shells. 

Key words. Stars: AGB and post-AGB - Stars: circumstellar matter - Stars: variable: general - Stars: evolution - Infrared: stars - 
Galaxy: bulge 


1. Introduction 

Low- and intermediate-mass stars pass along the asymptotic gi¬ 
ant branch (AGB) towards the end of their evolution. Stellar 
evolution models predict that stars approaching the tip of the 
AGB experience thermal pulses. In this so-called ‘thermal pulses 
AGB’ (TP-AGB) phase, stars steadily increase in luminosities 
and mass-loss rates. This increase is interrupted on timescales 
of several ten thousand years by short-term variations due to 
the pulses. The heavy mass loss leads to the formation of cir¬ 
cumstellar envelopes (CSE) of gas and dust. If the mass-loss 
rate surpasses M> 10~ 6 M o yr~\ the d ust s hell eventually be¬ 
comes opaque to visible light dHabinglll996h . In addition, the 
stars may change their photospheric chemistry from O- to C-rich 
as a consequence of the so-called third dredge-up following a 
thermal pulse. The range of initial masses, in which the star will 
end as C-rich, depends on the metallicity. This range, defined 
on th e main seq uence, is 2 $ Mms ~ 4 M co for solar metallicities 
dMarigo & Girardill2007i) . At lower metallicities, as for example 
in the Magellanic Clouds, this mass range is extended to lower 
masses, significant increasing the fraction of C stars in the AGB 
population at the end of the TP-AGB phase ( Blum et al.ll2006h . 
The stars ending the TP-AGB O-rich must have developed either 
from lo w-mass star s with main-sequence masses Mms ~ 2 M e 
dBertelli et al.ll2008l) . for which the TP-AGB phase is terminated 


before the conversion is completed, or from intermediate-mass 
stars with Mms ~ 4 M 0 , in which hot bottom burning (HBB) is 
active at the base of the convective envelope, preventing the for¬ 
mation of carbon. 

Low-mass stars are expected to experience only moderate 
mass-loss rates during most of their AGB lifetime, leading to the 
formation of optically thin CSE. Low-mass AGB stars have been 
discovered predominantly in the optical and near-infrared wave¬ 
length range. Classical Mira variables belong to these kind of 
stars. AGB stars descending from progenitors with larger main- 
sequence masses are expected to be more luminous and to have 
higher mass-loss rates. They probably spend a larger fraction 
of their TP-AGB life in an optically obscured state dMouhcinel 
l2002i) . which constrains their detection to the infrared or radio 
wavelength range. OH/IR stars, originally discovered by their 
radio OH maser emission and later identified in the infrared are 
an example of these kind of stars. OH maser emission is however 
also widespread among Mira variables, so that the term ‘OH/IR 
star’ encompasses low- and intermediate-mass stars. 

If low-mass stars develop high mass-loss rates at the end 
of their evolution on the TP-AGB, they will form optically 
thick CSEs at lower mass-loss rates than intermediate-mass stars 
because they have more compact shells with smaller outflow 
velocities. Also, because of their frequency they should out¬ 
number the intermediate-mass stars even among the obscured 
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stars. Howev er, in the Galactic disk, IChen et aD (1200ill and 
iLikkell (Il989h found that the reddest (e.g. with highest mass-loss 
rates) OH/IR stars are dominated by progenies of main-sequence 
intermediate-mass stars. A dearth of low-mass stars among in¬ 
frared bright sources in clusters of the Magellanic Clouds was 
also noted by Ivan Loon et al.l (120051) . Therefore, low-mass stars 
may never evolve into the high mass-loss regime observed for 
the very red OH/IR stars, or it might be a very brief phase be¬ 
cause the envelope mass above the core is rapidly depleted. 

These stars should have re lativ ely low luminosities, L < 8000 
L 0 for Mms ~ 2 M 0 (iBertelli et al .1120081) . Because of the un¬ 
certain distances in the Galactic disk and the strong variabil¬ 
ity, mean luminosities for AGB stars are difficult to determine 
and objects with relatively low luminosities are difficult to sin¬ 
gle out. A way to overcome this pitfall is to study AGB stars 
in the Ga lactic bulge, where a common distance can be used 
(Ivan der Veen & Habingfl990l) . AGB stars are luminous enough 
to be detectable there, especially in the mid-infrared. During 
their heavily obscured stage, they emit most of their radiation 
in the 10-30/rm range, which is less affected by the significant 
interstellar extinction towards many parts of the Galactic bulge. 

The Galactic bulge is the central part of the Milky Way 
galaxy, c onsisting of a dominant old population o f stars 
(t > 8 Gyr; IZoccali et al .1120031 : IVanhollebeke et ai1l2009h and a 
smaller population of younger s tars with no well constrained 
ages o f ~ 200 My r -_7 Gvr Avan Loon et ai]l2003l) or 1 - 3 Gyr 
dGroenewegen & Blommaerfl 20051) . to which the AGB stars in 
the Galactic bulge belong. In Baade’s window two metallic - 
ity populations have been found, of which the metal-poor pop¬ 
ulation is associated with the old stellar population and the 
metal-rich population is consistent with a younger stellar pop¬ 
ulation supporting a bar (ilHill et alJl20nl) . A met allicity gradi¬ 
ent was noted along the minor axis of the bulge ( IZoccali et all 
120081) with the metallicity decreasing towards higher latitudes. 
No lu minous C-rich AGB star s are known in the Galactic 
bulge felanco & Terndrutjll98 9). suggesting that the younger 
stellar population in the Galactic bulge may be restricted to 
low-mass stars with AGB evolutionary timescales too short to 
experie nce enough dredge-ups to convert them into C stars 
dBuelll 120131) . Then, the Galactic bulge OH/IR stars should 
also exclusively be descendants from low-mass stars, includ¬ 
ing the extremely red stars. This was implied by the findings 
of I van der Veen & Habingl (11990h . who derived low progenitor 
masses of 1.0 -1.4 M Q for a sample of these red sources from 
the IRAS survey. 

Although subsequent studies of AGB stars in the Galactic 
bulge also concluded that their samples originate from low- 
mass stars, luminosity distributions are not necessarily compat¬ 
ible with an ex clusive origin from low-mass stars. For exam¬ 
ple, jGTass|etl2 (2Q0l]) studied the variability properties of long- 
period, large-amplitude variable stars near the Galactic centre, 
and found a period-luminosity (P-L) relation that is compatible 
with the existence of a young population with initial masses in 
the range 2.5 -3.0 M e . In a later study on Mira variables in the 
optic al g ravitati onal lensing exper iment (OGLE) bulge fields, 
iGroenewegen & Blommaert (120051) found a different P-L rela¬ 
tion, which they explained with a population with low initial 
masses (1.5-2.0 M e ). However, possibly more massive stars 
with P>700 days, such as obscured OH/IR stars, mi ght have 
been m issed in the I-band of the OGLE survey. IBlommaert et al] 
( 200 61) studied a small sample of Galactic bulge AGB stars with 
near-infrared spectroscopy, and found a luminosity range from 
1700 to 7700 L 0 with a mean of 4100 L 0 , as expected for a pop¬ 
ulation of AGB stars with low initial masses. Again, their study 


was restricted to objects with relatively low mass-loss rates, and 
may therefore have excluded more massive stars. This can be in¬ 
ferred from lOiha et al.l (120071) . who focussed their study on high 
mass-loss AGB stars, and found a wider luminosity distribution 
with a peak at 8000 L 0 and an appreciable number of sources 
with a luminosity exceeding 20,000 L 0 . 

A numbe r of studies have been speci fically focus s ed on 
OH/IR stars. IBlommaert etafl (119981) and IWood et al.l (fk998l) 
found that the OH/IR stars o f the Galactic centre do not fit the 
classi cal Mira P-L relations dFeast et alJI 19891 : IWhitelock et al.l 
1199ll) . They systematically have lower luminosity and higher 
expansion velocity for a given period. OH/IR stars with periods 
as long as 1000 days, expected to have higher initial masses, 
were found to have low luminosities corresponding to Miras 
with masse s below 1.5 M 0 . These findings were explained by 
IWood et al .' ( 19981) with a higher metallicity in the Galactic cen¬ 
tre compared to the rest of the bulge. This fits into the general 
view that at low latitudes the metal-rich stellar population of the 
disk merges with the generally more metal-poor population of 
the bulge prevailing at larger latitude^ The origin from low-mass 
stars was confirmed bv lOrtiz et al.l (120021) for other OH/IR stars 
observed by the ISOGAL survey. They concluded that the small 
number of objects more luminous than M boi = -5.5 (L = 12,000 
L 0 ) indicates that the initial masses of stars in the sample are 
rarely higher than ~ 3 M 0 , but instead the masses are about 1-2 
M 0 , on average. 

All AGB star samples studied in the Galactic bulge were 
dominated by objects with blue or moderate reddened colours. 
This may have biased the studies towards low-mass objects, 
which would explain the scarcity of high-luminosity sources. 
In this paper, we study a sample of AGB stars of the Galactic 
bulge with extremely reddened colours. These stars are expected 
to be at the very end of the TP-AGB evolution. If low-mass stars 
evolve into the high mass-loss regime, the sample should include 
objects with rather low luminosities. If not, they must have orig¬ 
inated from intermediate mass stars with correspondingly higher 
luminosities. 

Currently, with new infrared surveys available, the spec¬ 
tral energy distributions (SED) of heavily obscured objects are 
covered with unprecedented coverage. Together with theoretical 
models of the CSEs of AGB stars, reliable bolometric fluxes can 
be derived without invoking uncertain bolometric corrections. 
Finally, it is also possible to obtain individual extinctions for 
each source leading to more reliable luminosities. 

The sample of extremely reddened AGB stars is presented in 
Sect.[2] SectionQ] describes the construction of the spectral en¬ 
ergy distributions (SED), their fitting with DUSTY models, and 
the estimation of luminosities and mass-loss rates. The nature of 
the sources in the sample is discussed in Sect.[4] and the conclu¬ 
sions are presented in Sect.0 


2. The sample 

Our sample of AGB stars of the Galactic bulge (hereafter: 
GB sampl e) was s elected from the sample of sources in 
iJimenez-Esteban et al.l (120061) . They studied around 100 IRAS 
sources distributed throughout the Milky Way with very red 
IRAS colours ([12] - [25] > 0.75 mag), located in a character¬ 
istic region of the IRAS two-colour diagram along the sequence 
of colo urs predicted for O-rich AGB stars with increasing mass 
loss dBediinl 11987h . and having a high (VAR >50) IRAS vari¬ 
ability index dBeichman et al.ll 1988h . The high value of this in¬ 
dex implies that all sources should be large-amplitude variables. 
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By comparison with the colo urs o f the ‘Arecibo sample’ of O- 
rich AGB stars ilimenez-Esteban et al.l 12005 1. they concluded 
that their sample of extremely reddened IRAS sources was also 
mainly made up of O-rich AGB stars highly obscured by their 
CSEs. To define the GB sample presented here, we selected the 
extremely reddened IRAS sources from iJimenez-Esteban et al] 
( 2006 1 with angular distances < 10° from the Galactic centre. 

Of the 39 sources selected in this region of the sky, two 
(IRAS 17317-3331 and IRAS 17411-3154) do not belong to 
the bulge. They are known to be foreground sources, located 
at a distance of 3.2 an d 1.2 kpc, respectively (ICohen et al.l 
l2005t [Yuasa et al.lll999l) . The remaining 37 sources, presented 
in TableQ] define our sample of extremely reddened AGB stars 
in the Galactic bulge. 

Although the sample is small, it is probably rather com¬ 
plete within its selection boundaries. Its flux density distribu¬ 
tion at 25 gm peaks at F 25 ~ 7 Jy, while the weakest source has 
f ’25 ~ 4 Jy (TableQ. AGB stars in the bulge matching the colour 
selection criterion [12] -[25] > 0.75 mag and F 25 > 7 Jy are prob¬ 
ably all contained in the sample, except those that are missing 
because of incompleteness of the IRAS catalogue. Surpassing 
the F 25 flux limit requires a minimum luminosity, which depends 
on the optical thickness of the CSE. This minimum luminosity 
is the lowest for sources in which the bulk of the emission is 
emitted in the 15-35 gm wavelength range. This minimum lu¬ 
minosity is also not biased by extinction, which is low at these 
wavelengths. As we show below, the minimum luminosity re¬ 
quired to surpass the F 25 flux limit is low enough for high opti¬ 
cal thickness, so that low-mass stars will also be present in the 
sample, if they develop the corresponding high mass-loss rates. 

In accordance with the two foreground disk sources identi¬ 
fied, we also expect contamination of the sample by background 
disk sources. The number of background disk sources can be 
roughly estimated as follows. Considering the foreground disk 
sources as “typical” for contaminating disk sources and using 4 
Jy as the F 25 flux limit, all such disk sources within the solar 
orbit around the centre of the Milky Way should have been de¬ 
tected by IRAS. Assuming the same disk source density in front 
and behind the Galactic bulge, the number of foreground and 
background disk sources is coarsely proportional to the surface 
area of the Galactic disk segments in front and behind the bulge. 
Using a distance to the Galactic centre of 8.0 (±0.5) kpc ( IReidl 
119931) . the adopted 10° angular size of the Galactic bulge cor¬ 
responds to a radius of ~ 1.4 kpc. Based on this consideration, 
we expect ~4 times more background sources than foreground 
sources. This yields about 8 sources out of 37, which could be 
background sources and whose luminosities would be underes¬ 
timated in the following. 

We searched in the literature for information on the pre¬ 
dominant chemistry in the CSE of our objects. Six sources 
have IRA S Low Resolution Spectra (LRS) spectral classifica¬ 
tions bv iKwok et ahl d 19971) . IRAS 17292-2727, IRAS 17316- 
3523, IRAS 17367-2722, IRAS 17418-2713, and IRAS 17495- 
2534 were classified as LRS class A, meaning that their 
spectra present the 9.1 pm silicate dust feature in absorption. 
IRAS 18092-2347 was classified as LRS class I, meaning that 
the LRS was noisy or incomplete. For the rest of the sources, we 
searched for LRS spectra in the Calgary databas^J Only the LRS 
of IRAS 17367-3633 was found and visually classified. It shows 
a clear 9.1 pm silicate dust feature in absorption. We also visu¬ 
ally inspected the LRS of IRAS 18092-2347, which also_showed 
the silicate dust feature in absorption. In addition, iGolriz et~akl 


1 http://www.iras.ucalgary.ca/saldatabase/index.html 


(l2014lf have recently observed three of our sources that lack 
LRS classifications (IRAS 17251-2821, IRAS 17276-2846, 
and IRAS 17323-2424) with the Spitzer- Infrared Spectrograph. 
All three Spitzer spectra show the 9.1pm silicate in absorption. 
We classified these five sources as LRS class A. 

In addition, all sources of our sample except 
IRAS 17128-3528 were observed for OH maser emission 


( te Lintel Hekkert et al.ll 1 99 it iDavid et ail 1 993flSevenster et al.l 


11997 ). Of these, eight were not detected, two (IRAS 17350- 
2413 and IRAS 18019-3121) presented single-peak OH 
emission, and the rest showed the typical double-peak spec¬ 
trum commonly shown by OH/IR stars. Expansion velocities 
of the CSEs were derived for our objects from the velocity 
separations between the two OH peaks. For IRAS 17495-2534, 
with no OH maser emission detected, the expansion velocity 
was obtained from CO data dLoup et al.l 11993 ). Four sources 
(IRAS 17207-3632, IRAS 17392-3020, IRAS 17418-2713, 
and IRAS 18092-2347) were observed fo r H 2 O maser emission 
dDeacon et al.l 120071: ISuare zetal. 2007 ). On ly IRAS 17392- 
3020 was detected by Deacon et al j ( 20071) in April 2003, 
although it was not detected two years later, in September 2005, 
by [Suarez et al.l (120071) . Finally, two sources (IRAS 17292- 
2727 and IRAS 1 7341-3529) were detected a s SiO masers 
(iNvman et al.lll998i : lNakashima & Deguchi82003l) . Evidence for 
the O-rich nature of 30 sources is therefore provided by either 
detection of a maser and/or the presence of the silicate dust 
feature. 

In TableQ] we present the observational parameters for the 
sources in our sample and summarise the above information. 
The IRAS object name is listed in the first column. Columns 
2 to 5 give the equatorial and Galactic WISE coordinates (see 
Sect l3.fl ). Columns 6 and 7 give the IRAS colours. The IRAS 
LRS spectral classification and the IRAS variability index are 
given in the eight and ninth columns, respectively. In columns 10 
to 12, we summarise the detection status for OH, HiO, and SiO 
maser emission, respectively. The expansion velocity, if avail¬ 
able, is listed in column 13, and the estimated individual ex¬ 
tinction in the line of sight used to de-redden the observational 
photometry (see Sect. 13. 21 ) is in column 14. Some comments on 
individual SEDs, discussed below, are given in the last column. 


3. Analysis 

3. 1. Spectral energy distributions 

To create the SEDs of the sources of our sample we col¬ 
lected, when available, photometric information from astronom¬ 
ical catalogues. To do that, we took advantage of the Virtual 
ObservatonQ] (VO). In particular, we made use of the Multiple 
Cone Search utility of Topca0, an interactive graphical viewer 
and editor for tabular data that allows the user to examine, anal¬ 
yse, combine, and edit astronomical tables. With Topcat, we 
cross-matched our sample with several available astronomical 
catalogues within the VO, using different search radii depending 
on the particular catalogue. 

This way, we populated the SEDs from the near- to the far- 
infrared, wavelength range where these objects radiate the ma¬ 
jority of their energy. We collected data from the following cata¬ 
logues: Two Micron All Sky Survey (2MASS: ICutri et aljl2003l) 
and VISTA variables in the Via Lactea catalogue (VVV-DR1; 
ISaito et ahll2012l) at 1.25, 1.65, and 2.2pm (filter J, H, and K s , 

2 http://www.ivoa.net 

3 http://www.star.bris.ac.uk/~mbt/topcat/ 
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Table 1 . The GB sample. Observational parameters. 




(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

GO) 

(11) 

(12) 

(13) 

(14) 

(15) 

IRAS 


WISE coordinates 



IRAS 



Radio data' 


V expf 

&-Ks 

Comment® 

name 

RA(J2000) 

Dec(J2000) 

Gal Long (°) 

Gal Lat (°) 

[12]-[25f 

[25]-[60]“ 

LRS* 

Var' 7 

OH 

h 2 o 

SiO 

[km s _1 ] 

mag 


16582-3059 

17:01:29.85 

-31:04:19.1 

352.843477 

6.688608 

0.89 

-0.86 


57 

Y 



13.4 

0.15 


17030-3053 

17:06:14.06 

-30:57:38.2 

353.547542 

5.945143 

0.77 

-0.73 


99 

Y 



15.3 

0.15 


17107-3330 

17:14:05.05 

-33:33:57.9 

352.421186 

3.067547 

0.86 

-0.96 


99 

N 




0.40 


17128-3528 

17:16:12.91 

-35:32:12.8 

351.072561 

1.563949 

1.52 

0.12 


98 





0.82 

NIR-Exc 

17151-3642 

17:18:29.53 

-36:46:01.9 

350.333858 

0.477947 

0.99 

-0.30 


99 

N 




1.63 


17171-2955 

17:20:20.35 

-29:58:21.7 

356.131283 

4.050980 

0.82 

-1.30 


99 

N 




0.30 


17207-3632 

17:24:07.27 

-36:35:40.6 

351.118359 

-0.351526 

0.95 

0.46 


99 

Y 

N 


16.1 

2.45 


17251-2821 

17:28:18.61 

-28:24:00.5 

358.413518 

3.491327 

0.94 

-0.75 

A c 

99 

Y 



16.0 

0.35 


17276-2846 

17:30:48.32 

-28:49:02.0 

358.366433 

2.804559 

1.25 

-0.12 

A c 

79 

Y 



15.5 

0.48 

PPN C ; DP 

17292-2727 

17:32:23.57 

-27:30:01.2 

359.662488 

3.230668 

0.81 

-0.79 

A 

64 

Y 


Y 

17.6 

0.53 

FIR-Exc 

17316-3523 

17:34:57.50 

-35:25:52.5 

353.298313 

-1.536651 

1.02 

-0.72 

A 

99 

Y 



12.3 

0.98 


17323-2424 

17:35:25.92 

-24:26:30.5 

2.611416 

4.308932 

0.99 

-0.73 

A c 

97 

Y 



13.8 

0.60 


17341-3529 

17:37:30.29 

-35:31:04.4 

353.504370 

-2.020286 

0.88 

-1.26 


97 

N 


Y 


0.50 


17350-2413 

17:38:08.86 

-24:14:49.2 

3.107790 

3.890090 

0.85 

-1.50 


99 

Sp 




0.32 


17351-3429 

17:38:26.28 

-34:30:40.7 

354.457640 

-1.643595 

0.95 

-0.58 


99 

Y 



20.0 

0.56 


17361-2358 

17:39:14.82 

-23:59:56.8 

3.451800 

3.809114 

0.80 

-0.88 


99 

N 




0.31 


17367-2722 

17:39:52.43 

-27:23:31.7 

0.646949 

1.889337 

0.98 

-0.47 

A 

97 

Y 



7.7 

0.54 

FIR-Exc 

17367-3633 

17:40:07.62 

-36:34:41.4 

352.888394 

-3.033409 

0.88 

-0.75 

A c 

92 

Y 



14.5 

0.46 

FIR-Exc 

17368-3515 

17:40:12.98 

-35:16:40.8 

354.002057 

-2.360174 

0.77 

-0.71 


99 

Y 



13.5 

0.35 


17392-3020 

17:42:30.54 

-30:22:07.4 

358.425640 

-0.174669 

1.09 

0.81 


54 

Y 

Y,N 


21.1 

1.80 


17418-2713 

17:44:58.73 

-27:14:42.7 

1.369081 

1.002630 

1.34 

-0.12 

A 

99 

Y 

N 


15.2 

0.76 


17428-2438 

17:45:56.94 

-24:39:57.9 

3.685349 

2.159567 

0.86 

-1.04 


73 

Y 



13.8 

0.58 


17495-2534 

17:52:39.54 

-25:34:39.1 

3.684492 

0.388289 

0.80 

-1.31 

A 

99 

N 



16.07 

1.28 


17504-3312 

17:53:50.25 

-33:13:26.5 

357.219129 

-3.708828 

1.17 

0.21 


84 

Y 



15.7 

0.26 

PPN C ; DP 

17521-2938 

17:55:21.80 

-29:39:13.0 

0.471972 

-2.191825 

1.07 

-0.17 


98 

Y 



16.3 

0.28 


17545-3056 

17:57:48.41 

-30:56:25.8 

359.621118 

-3.292745 

1.36 

-0.04 


89 

Y 



15.0 

0.29 


17545-3317 

17:57:49.20 

-33:17:47.6 

357.573330 

-4.466135 

0.86 

-0.91 


99 

Y 



13.9 

0.25 


17583-3346 

18:01:39.29 

-33:46:00.5 

357.556949 

-5.393255 

0.76 

-0.99 


99 

Y 



13.4 

0.20 

FIR-Exc 

17584-3147 

18:01:42.04 

-31:47:54.7 

359.284968 

-4.439273 

1.04 

-0.53 


99 

Y 



16.8 

0.20 


18019-3121 

18:05:12.26 

-31:21:43.2 

0.032084 

^1.878452 

0.99 

-1.00 


99 

Sp 




0.21 


18040-2726 

18:07:08.92 

-27:25:52.4 

3.685892 

-3.345684 

1.05 

-0.58 


84 

Y 



16.3 

0.19 


18040-2953 

18:07:18.44 

-29:53:11.2 

1.547809 

-4.562185 

0.82 

-0.64 


99 

N 




0.17 


18091-2437 

18:12:16.15 

-24:36:42.9 

6.715421 

-2.997292 

0.80 

-0.60 


99 

Y 



15.7 

0.47 


18092-2347 

18:12:20.44 

-23:46:55.9 

7.452681 

-2.614618 

0.94 

-0.15 

I-A c 

92 

Y 

N 


17.4 

0.45 


18092-2508 

18:12:21.87 

-25:07:20.6 

6.276710 

-3.260295 

0.80 

-0.80 


91 

Y 



14.5 

0.40 

NIR-Exc 

18195-2804 

18:22:40.21 

-28:03:08.8 

4.760139 

-6.667509 

0.81 

-0.70 


99 

Y 



16.3 

0.13 


18201-2549 

18:23:12.28 

-25:47:58.7 

6.828182 

-5.738816 

0.77 

-1.50 


98 

N 




0.16 



Notes: " IRAS colours as defined in ljimenez-Esteban et al.l (1 20051) : b IRAS Low Resolution Spectra; c Classification by this work (see Sections l2land !3. lt : rf IRAS variability index. c Radio data. Y: 
detection, N: no detection, Sp: single peak; f Expansion velocity from the OH maser velocity range, or from CO data for IRAS 17495-2534; g SED peculiarities (see Section [3Tt . DP: double 
peaked, NIR-Exc: near-infrared excess, FIR-Exc: far-infrared excess, PPN: proto-planetary nebula candidate. 
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respectively); GLIMPSE ( Churchwell et al.ll2009l) at 3.56, 4.51, 
5.76, and 7 .96 um (filter s II, 12 ,13, and 14, respectively); WISE 
catalogue (Wright et al.l l2010h at 3.4, 4.6, 11.6, and 22.1yum 
(filters Wl, W2, W3, a nd W4, respective ly); MSX6C Infrared 
Point Source Catalogue dEgan et al.ll2003l) at 8.28, 12.13, 14.65, 
and 21.3 fim (filters A, C, D, and E, respectively); AKARI 
Point Source Catalogue (llshihara et al.l 120101 : lYamamura et alJ 
l20K)l) at 9, 18, 65, 90, 140, and 160^ (filters S09, S18, S65, 
S90, S140, and S160, resp ectively); and IRAS Point Sources 
Catalogue dBeichman et al.lU988l) at 12, 25, 60 and 100 /jm (fil¬ 
ters 12,25,60, and 100, respectively). Since WISE data were not 
available through the VO, they were downloaded directly from 
the NASA/IPAC Infrared Science Archive. 

We plotted the SED for each source and visually in¬ 
spected them to check the consistency of the photometric data. 
When d eviat ing photometric data were found, we used AladirQ 
dBonnarel et al.1 12000 1. another VO tool, to visualise the astro¬ 
nomical images of the field and to check the reliability of the 
selected counterpart. Aladin allows users to visualise and anal¬ 
yse digitised astronomical images, and superimpose entries from 
astronomical catalogues or databases available from the VO ser¬ 
vices. 

We found few inconsistencies, basically because of chance 
coincidences with field stars within our search radii. In all these 
cases the offending counterparts were removed. In the case of 
missing counterparts both in 2MASS and VISTA catalogues, we 
used a nearby (<30") faint star to define an upper flux limit, 
preferably in the VISTA catalogue. 

The photometric data are presented in Table[2] and the SEDs 
are plotted in Fig. 7. This figure shows the observational data 
after correcting for interstellar extinction and the DUSTY model 
that best fits the observations (see Sect. l3.2l >. Most of the sources 
present the typical SEDs of an AGB star heavily obscured by an 
optically thick CSE, with the maximum of their emission in the 
mid-IR. Note that most of the sources were not detected in the 
near-IR, or only barely in the K, band. 

Two sources (IRAS 17276-2846 and IRAS 17504-3312) 
show a peculiar double-peaked SED (Fig. 7). In the case of IRAS 
17276-2846, the photospheric emission extends into the visual 
so that we added VVVDR1 data for filters Z (0.88 yum) and Y 
(1.02 /rm) to the SED. These SEDs appear as soon as the mass- 
loss rates have greatly decreased and the obscuration has sub¬ 
stantially diminished. The bluer peak at ~ 2 //m is then attributed 
to the reddened stellar photosphere, and the redder peak corre¬ 
sponds to the emission at longer wavelengt hs by the cool cir- 
cumstellar dust (e.g. see Fig. 9 in lSanchez Contreras et al.l2008l) . 
This kind of drop in mass-loss rate can be expected in the after- 
math of a thermal pulse or when the stellar envelope has nearly 
been lost at the end of AGB evolution. Stars with this character¬ 
istic SED are generally considered to be proto-planetary nebula 
(PPN) dOudmaiierl 119961) . and therefore we assume that these 
stars have left the AGB. Since post-AGB evolution is outside 
the scope of the present paper, the two PPN candidates are not 
included in the following analyses. 

3.2. Observed bolometric flux 

The bulge sample of AGB stars is nowadays covered over a large 
part of their SED by a multitude of infrared surveys. This allows 
us to determine the bolometric flux with much higher accuracy 
compared to earlier times, when (uncertain) bolometric correc¬ 
tions or (unsuitable) black-body fitting had to be applied to cover 

4 http://aladin.u-strasbg.fr/ 


the unobserved parts of the SED (e.g. Ivan der Veen & Breukersl 
1 989 : Ortiz et all 12002 4 However, the photometric data avail¬ 
able is not complete for all sources so that a determination of 
the bolometric flux by integration of an interpolated curve be¬ 
tween the available data points will be unsatisfactory. Also bolo¬ 
metric corrections (although small ones) for the spectral ranges 
outside the infrared are still required. We therefore decided to fit 
the observed SEDs with models of stars obscured by a dusty ra¬ 
dial symmetric CSE and to determine the bolometric fluxes from 
the models. Because of the large distance of the Galactic bulge, 
the stars may be observed through large columns of intervening 
dust, such that extinction may have a noticeable influence on the 
shape of the SEDs and on the observed bolometric flux. 

To de-redden the photometric data, we corrected the flux at 
2.2 /urn (Ks band) using the extinction maps towards the Galactic 
bulge given by iGonzalez et al.l (120121) . The absolute extinc¬ 
tions towards the sources of our sample are 0.13 < A^ < 2.45 
mag (see Table Q}. Two sources (IRAS 16582-3059 and 
IRAS 17030-3053) are not covered by the maps and we assigned 
Aks = 0.15 mag to them. About half of th e sou rces are also cov¬ 
ered by the extinction maps of iNidever et al. ( 20121) . Their ab¬ 
solute extinctions are 0.14 + 0.04 mag higher for A^ s <1.0 mag 
than those of IGonzalez et al.l (120121) . Because the uncertainties 
given in both papers are 0.1 mag, we considered this difference 
as not significant. For higher extinctions, the scatter <x between 
the maps is larger with cr = 0.3 mag. 

To correct the other photometric bands, we determined the 
ex tinction coefficient s Aj/A^., from the extinc tion curve s given 
bv lGhen et al.l (12013l) up to A — 8 yum and from lGao et alJ (120091) 
up to 24 /urn. The extinction coefficients beyond A = 24 fjm were 
determined by extrapolation using a power law with yd = -1.7. 

The de-reddened SEDs were jnodelled using the radiation 
transport code DUSTY (Ivezic et al.l 1 19991) . Because the ob¬ 
served photometry taken at different epochs is likely affected by 
variability, the photometric data are expected to scatter around 
the mean SED. The relatively sparse coverage of the SEDs with 
photometry and the variability preclude the determination of ac¬ 
curate models. We therefore decided to use a standard model, 
in which most of the parameters are fixed, and only the optical 
depth of the shell was left as variable parameter. 

Combining the information from the IRAS LRS classifi¬ 
cations and the maser detections, O-rich chemistry can be in¬ 
ferred for 28 o bjects. Carbon stars are e xtremely rare in the 
Galactic bulge (Blanco & Terndrup 1989 ). Hence, we decided 
to assume O-rich chemistry for the remaining seven objects, 
following iLewisI (119921) . who found that infrared sources with¬ 
out OH maser detection but with IRAS colours similar to those 
of OH/IR stars are most probably O-rich variable AGB stars. 
Consequently, we only used models with O-rich chemistry. 

The standard model assumes a central star with an effec¬ 
tive temperature T e yy = 2500 K, and a dust condensation tem¬ 
perature T d = 1000 K. We use d the optic al constants for amor¬ 
phous cold silic ates from Ossenko pf et al.l (1992 ) and took the 
standard MRN (Mathis et al.l 119771) dust size distribution with 
n{a) oca~ 3 ' 5 , where n is the number density and a is the size of 
the grains. The grain sizes were limited to 0.005 < a < 0.25 yum. 
The inner radius of the CSE is determined by DUSTY from 
the dust condensation temperature assumed. The outer radius of 
the shell was set to R max = 100 Instead of describing the 
density distribution in the shell by an analytic profile, the option 
to calculate the wind structure from hydrodynamics was cho¬ 
sen, which_allows_us todetermine mass-loss rates from the mod¬ 
els (Ivezic & Elitzuriri995l) . Models were calculated for optical 
depths 0.2 < Tg j^m < 50 and the best model was found by min- 
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Table 2. Photometric data collected. 


ON 


Filter 

J 

H 

K s 

W1 

11 

12 

W2 

13 

14 

A 

S09 

W3 

12 

C 

,1 (pm) 

1.25 

1.65 

2.20 

3.4 

3.56 

4.51 

4.6 

5.76 

7.96 

8.28 

8.61 

11.6 

11.60 

12.13 

Flux in Jy x 

10~ 4 

ur 4 

10 -4 

10" 3 

10“ 3 

10“ 3 

10“ 2 

10^ 2 

lO^ 2 

io-' 

10-' 

io-' 

10° 

10° 

16582-3059 

"2.6 

a 3.9 

“25.1±1.3 

138±3 



77.3±1.6 




20.5±0.2 

23.1 ±0.4 

2.5±0.2 


17030-3053 

"1.6 

“1.8 

“66.8±1.7 

514±16 



310±10 



20.0±0.8 

31.3± 1.6 

32.0±0.7 

2.0±0.3 

3.2 

17107-3330 

"4.0 

“29.3±1.5 

“609±16 

680±20 



270±9 



20.3±0.8 

47±3 

34.2±0.7 

2.60±0.16 

2.87±0.16 

17128-3528 

“24 

*1.52±0.14 

*17.6±0.2 

25.0±1.0 



20.9±0.4 



29.0±1.2 

50±40 

36.3±0.7 

2.2±0.3 

7.9±0.4 

17151-3642 

*0.2 

*3.5 

*0.7 

132±4 



69.0±1.6 



47±2 

20.3±0.2 

30.4±0.6 

5 ±2 

7.1 ±0.4 

17171-2955 

*0.5 

*20.2 

*156.4±0.3 

219±5 



100±3 



35.2±1.4 

23.0±0.4 

17.5±0.3 

3.2±1.7 

4.7±0.3 

17207-3632 

*0.2 

*1.39 

*38.9 

43.7±1.6 



163 ±5 



56±2 

54.7±0.7 

104.7±1.2 

6.2± 1.2 

8.8±0.4 

17251-2821 

"75 

*12.1 

*1.2 

38.0±1.4 

192±6 

859±40 

41.8±0.8 

263±4 

404±9 

21.3±0.9 


24.2±0.5 

3.6±0.8 

3.7±0.2 

17276-2846'' 

“700±30 

“1700 

*1441.0±1.3 

135±4 

138±4 

261 ±7 

32.7±0.7 

137±3 

340±20 

21.3±0.9 


34.0±0.7 

2.5±0.5 

3.4±0.2 

17292-2727 

“70 

“64 

“63 

870±30 



780±30 



80±3 

88±12 

145.1±0.8 

11.6±0.6 

11.4±0.6 

17316-3523 

*0.7 

*1.1 

“42±2 

*24.4±0.3 

437±12 



608±16 



62±3 

61±7 

145.1±0.5 

10±3 

10.7±0.5 

17323-2424 

*0.4 

*3.8 

*0.6 

24.2±0.9 



44.1 ±0.9 



23.4±1.0 


30.6±0.6 

3.4±0.5 

2.9±0.2 

17341-3529 

“22.6±1.3 

“316±8 

“1590±30 

1140±40 



409±14 



83±3 

41±4 

63.7±0.8 

6.3±1.6 

12.5±0.6 


*3.0±0.3 

*67.0±0.4 

*279.2±0.5 












17350-2413 

*0.7 

*0.8 

*175.7±0.3 

178±4 



63.8± 1.3 



14.8±0.6 


12.9±0.2 

3.6±1.3 

1.99±0.13 

17351-3429 

*1.3 

*1.8 

“107±4 

*16.3±0.5 

414± 11 



457±15 



45 ±2 


58.5±0.7 

3.2±0.7 

6.1 ±0.3 

17361-2358 

*0.5 

*0.8 

*19.2±0.3 

253±7 



253±7 



19.6±0.8 

63.5±0.8 

49.2±0.4 

3.6±0.8 

2.70±0.15 

17367-2722 

*1.0 

*1.1 

*1.1 

29±2 



50±1 



24.1 ±1.0 


32.0±0.7 

3.7±0.4 

3.8±0.2 

17367-3633 

*0.3 

*8.4±0.2 

“66±2 

*43.9±0.4 

690±20 



710±20 



62±3 


147.2±1.0 

8.8±1.1 

8.2±0.4 

17368-3515 

*0.2 

*1.2 

*0.8 

27.9±1.6 



60.6±1.3 



31.4± 1.3 


39.3±0.7 

4.0± 1.1 

4.6±0.2 

17392-3020 

*0.2 

*0.3 

*1.0 

74±5 

37.7± 1.5 

350±20 

80.7±1.8 

179±3 


40.1 ±1.6 

20.30±0.03 

33.9±0.7 

6.3±0.5 

5.2±0.3 

17418-2713 

*1.0 

*1.3 

*2.0 

40±3 

22.6±1.0 

490±30 

175±4 

368±7 


179±7 

179.3±0.6 

384.4±1.1 

15±3 

29.8± 1.5 

17428-2438 

*0.9 

*1.1 

*20.6±0.4 

292±8 

246±8 

630±40 

308±9 

240±5 


36.6± 1.5 

45± 15 

57.0±0.6 

6.6±0.4 

4.73±0.13 

17495-2534 

*0.3 

*0.3 

*144.3±0.5 

1108±49 

370±20 

8400±800 

1140±30 

436±11 


85 ±4 


314.5± 1.2 

29±5 

13.7±0.7 

17504-3312 

*54.1±0.3 

“70.3±1.4 

“209±4 

59±2 



11.4±0.3 



7.3±0.3 

6.97±0.05 

18.1 ±0.3 

4.1±0.8 

2.42±0.15 



* 193.1 ±0.4 

*425.8±0.4 












17521-2938 

“31 

“44 

*2.3 

3.6±1.7 

29.2±1.5 

269±9 

16.5±0.4 

130±3 


15.2±0.6 

26±6 

18.3±0.4 

3.1±0.6 

2.36±0.12 

17545-3056 

*1.3 

*1.4 

*1.6 

2.2 

1.52±0.08 

31.6± 1.2 

5.72±0.16 

33.4±0.7 

119±3 

11.1 ±0.5 


22.4±0.4 

2.9±0.4 

2.53±0.13 

17545-3317 

*0.8 

*1.3 

“45±3 
*11.4±0.3 

204±5 



140±4 



46±2 


32.7±0.6 

3.6±0.3 

6.2±0.3 

17583-3346 

*25.8±0.2 

“39±2 

“522±9 

631±19 



210±6 





36.2±0.7 

2.00±0.14 




*265.1±0.2 

*1082.0±1.0 












17584-3147 

*0.7 

*0.7 

*1.0 

46.9±1.6 



95 ±2 



22.6±0.9 


67.7±0.8 

6.9±1.0 

3.4±0.2 

18019-3121 

*0.5 

*0.7 

*0.8 

25.2±0.9 



50.0±1.0 



24.5±1.0 

28±4 

29.3±0.6 

3.5±1.1 

3.0±0.2 

18040-2726 

*1.1 

*1.1 

*0.8 

39.1± 1.3 



97±2 



24.2±1.0 

51 ± 15 

59.4±0.9 

6.5±0.5 

3.5±0.2 

18040-2953 

*0.9 

*1.0 

*8.6±0.3 

345±9 



492±17 



28.3± 1.2 

30±4 

84.4±0.6 

5.0±0.8 

3.7±0.2 

18091-2437 

*0.8 

*0.7 

*10.6±0.3 

177±4 

236±8 


150±4 

212±5 


16.3±0.7 

20.2±0.9 

32.8±0.7 

4.6±0.7 

2.30±0.14 

18092-2347 

*0.5 

*1.0 

*1.0 

4.3± 1.1 

13.9±0.6 

279±13 

36.5±0.7 

230±5 


50±2 

55±3 

121.2±0.9 

10.9±0.8 

9.8±0.5 

18092-2508 

*1.5±0.2 

*2.4±0.3 

*20.2±0.3 

311±9 



352±13 



49±2 

51.8±7 

58.4±0.5 

5.6±0.7 

6.6±0.3 

18195-2804 

“0.4 

*0.5 

*21.76±0.16 

362±8 



341±10 




20.6±6 

49.2±0.5 

3.6±0.8 


18201-2549 

“11.1 ±0.7 

“323±7 

“3260±60 

2340±140 



678±16 




93.627±0.003 

95.0±0.5 

5.0±0.7 



*14.22±0.18 

*366.9±0.3 

*2759±3 













Notes: Data without errors are upper limits; " 2MASS data; * VVV-DR1 data ; c No confirmed detection; d VVV-DR1 fluxes in Jy x 10 5 : F ; =761±2 and Fy=3017±4. 
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Table 2. Photometric data collected. Continued 


Filter (A ( pm )) 

A ( pm ) 

Flux in Jy x 

D 

14.65 

10 ° 

S 18 

18.39 

io-' 

E 

21.3 

10 ° 

W 4 

22.1 

io-' 

25 

23.88 

10 ° 

60 

60 

10 ° 

S 65 

65 

10 ° 

S 90 

90 

io- 1 

100 

100 

10 ° 

SI 40 

140 

10 ° 

S 160 

160 

10 ° 

16582-3059 




51 . 7 ± 0.9 

5 . 8 ± 0.6 

2 . 6 ± 0.3 

c 1.6 

13±2 

20 


c 1.3 

17030-3053 

3 . 5 ± 0.2 


5.5 

58 . 0 ± 0.8 

4.1 ± 0.7 

2 . 1 ± 0.2 



26 



17107-3330 

3 . 4 ± 0.2 


3 . 8 ± 0.2 

64 . 9 ± 0.8 

5 . 7 ± 1.0 

2 . 4 ± 0.5 



20 



17128-3528 

13 . 0 ± 0.8 

93 . 0 ± 0.3 

13 . 8 ± 0.8 

114 . 9 ± 1.5 

9 . 1 ± 1.4 

10 . 1 ± 0.8 

17±5 

86±8 

55 

c 9.9 

c 7.3 

17151-3642 

11 . 2 ± 0.7 


11 . 3 ± 0.7 

73 . 2 ± 1.6 

13±4 

10.1 ± 1.6 

c 6.8 

71 ±5 

320 

c 30±4 

c 82±7 

17171-2955 

5 . 8 ± 0.4 


6 . 1 ± 0.4 

33 . 6 ± 0.3 

7±4 

2 . 1 ± 0.6 

c l.l 

16 . 5 ± 0.2 

18 



17207-3632 

15 . 4 ± 0.9 


12 . 6 ± 0.8 

194±3 

15±3 

23±3 

16±4 

113±6 

1500 


c 45±9 

17251-2821 

4 . 9 ± 0.3 

47±8 

5 . 8 ± 0.4 

58 . 8 ± 0.7 

8 . 5 ± 1.6 

4 . 2 ± 0.4 

c 1 . 7 ± 0.6 

22±2 

20 

c 1.8 


17276 - 2846 rf 

5 . 7 ± 0.4 

59 . 20±4 

5 . 4 ± 0.3 

75 . 7 ± 0.9 

7 . 9 ± 0.9 

7 . 2 ± 0.8 

c 6.2 

37±3 

30 

c 1.4 


17292-2727 

16.1 ± 1.0 

153 . 3 ± 0.7 

19 . 5 ± 1.2 

213±6 

24±2 

11 . 8 ± 1.1 

8 . 9 ± 0.4 

71 ±9 

12 

4 . 6 ± 1.2 

c 4.4 

17316-3523 

16 . 4 ± 1.0 


19 . 6 ± 1.2 

246±4 

25±4 

13±2 



170 



17323-2424 

5 . 5 ± 0.3 

43 . 0 ± 1.6 

7 . 1 ± 0.5 

67 . 9 ± 0.9 

8 . 5 ± 1.7 

4 . 3 ± 0.5 

c 3.8 

19±4 

13 


c 0.66 

17341-3529 

13 . 3 ± 0.8 

87±10 

16 . 4 ± 1.0 

116 . 9 ± 1.3 

14±2 

4 . 5 ± 0.9 



100 



17350-2413 

2 . 76 ± 0.18 


3.1 

27 . 6 ± 0.3 

7 . 8 ± 0.6 

2 . 0 ± 0.2 

c 1.4 

13.1 ± 0.6 

13 

c 0.61 


17351-3429 

8 . 6 ± 0.5 

106 . 9 ± 1.6 

10 . 2 ± 0.6 

117.1 ± 1 . 1 

7 . 8 ± 1.4 

4 . 6 ± 0.7 

6 . 0 ± 0.8 

29±6 

160 

c 2.0 

c 0.13 

17361-2358 

3 . 7 ± 0.2 


4 . 1 ± 0.3 

85 . 9 ± 1.7 

7 . 5 ± 1.4 

3 . 3 ± 0.4 



11 



17367-2722 

6.1 ± 0.4 

90±30 

6 . 6 + 0.4 

76 . 7 ± 1.3 

9 . 2 ± 1.2 

6 . 0 ± 0.5 

c 4.0 

44±7 

36 

c 4 . 4 ± l.l 

c 4.1 

17367-3633 

11 . 6 ± 0.7 


13 . 3 ± 0.8 

204±4 

20±2 

9 . 9 ± 1.0 

11 ±5 

69 . 6 ± 0.4 

10 

c 6±4 

c 0.92 

17368-3515 

7 . 3 ± 0.5 

74 . 3 ± 1.3 

7 . 5 ± 0.5 

85 . 0 ± 1.0 

8±2 

4 . 2 ± 1.0 

c 5 . 9 ± 1.5 

21 ±3 

13 



17392-3020 

9 . 0 ± 0.5 

90±30 

9 . 4 ± 0.6 

91 . 6 ± 1.6 

17 . 3 ± 1.6 

36±7 



210 



17418-2713 

58±4 


68±4 

720±9 

51 ±7 

46±5 



500 



17428-2438 

6 . 3 ± 0.4 

103 . 5 ± 1.1 

7 . 2 ± 0.4 

110 . 5 ± 1.7 

14 . 6 ± 1.2 

5 . 6 ± 0.6 



47 



17495-2534 

19 . 9 ± 1.2 

500±40 

23 . 7 ± 1.4 

521±8 

61±8 

18±2 

17±4 

88 . 9 ± 1.7 

370 



17504-3312 

4 . 7 ± 0.3 

59 . 34 ± 0.13 

7 . 8 ± 0.5 

87 . 6 ± 1.0 

12 . 0 ± 1.2 

15±2 



43 



17521-2938 

4 . 3 ± 0.3 

69±13 

5 . 4 ± 0.3 

49 . 8 ± 0.7 

8 . 2 ± 1.7 

7 . 0 ± 1 . 1 



110 



17545-3056 

5 . 5 ± 0.3 

50±20 

6 . 5 ± 0.4 

75 . 6 ± 0.9 

10 . 2 ± 0.9 

9 . 8 ± 1.2 

5 . 3 ± 1.2 

c 45.56 

57 

c 7.0 

c 0.25 

17545-3317 

7 . 9 ± 0.5 


9 . 1 ± 0.6 

50.1 ± 0.5 

8 . 0 ± 1 . 1 

3 . 5 ± 0.4 

c 1.7 

14±2 

34 


c 0.42 

17583-3346 


35 . 0 ± 0.8 


59 . 0 ± 0.6 

4 . 0 ± 0.3 

1 . 6 ± 0.4 

c l.l 

13 . 0 ± 1.6 

26 

c 1.0 


17584-3147 

5 . 1 ± 0.3 

99±17 

6 . 6 ± 0.4 

167 ±4 

18±3 

11 . 0 ± 1.3 

6 . 82 ± 0.13 

c 43.67 

12 



18019-3121 

5 . 7 ± 0.4 

69 . 17 ± 0.04 

5 . 2 ± 0.3 

64 . 5 ± 0.8 

9±2 

3 . 4 ± 0.9 



62 



18040-2726 

6 . 5 ± 0.4 

130±20 

6 . 4 ± 0.4 

137±2 

17 . 0 ± 1.4 

9 . 9 ± 1.2 

8 . 8 ± 0.4 

69±3 

66 

c 1.7 

c l.l 

18040-2953 

5 . 7 ± 0.4 

60±7 

6 . 3 ± 0.4 

143±2 

11±2 

5 . 9 ± 1.0 

c 2.9 

25 . 01 ± 0.02 

42 


c 1.7 

18091-2437 

3 . 1 ± 0.2 


3 . 8 ± 0.2 

63 . 8 ± 0.7 

9 . 7 ± 1.5 

5 . 6 ± 0.6 



86 



18092-2347 

17 . 7 ± 1.1 


1 8 . 1 ± 1 . 1 

245±5 

26±2 

22 . 4±3 

13 . 2 ± 0.3 

115±7 

110 


C 4J 

18092-2508 

9 . 5 ± 0.6 

121 . 0 ± 0.9 

10 . 6 ± 0.7 

115±2 

11 . 8 ± 1.1 

5 . 6 ± 0.7 

C 4J 

26 . 9 ± 0.9 

90 

c 0.79 

c 2.0 

18195-2804 


30 . 4 ± 0.5 


95.1 ± 1.7 

7 . 7 ± 1.5 

4 . 0 ± 0.5 

1 . 7 ± 0.3 

16±2 

17 

c 0 . 6 ± 0.8 

c 1.4 

18201-2549 


64 . 9 ± 0.4 


133 . 7 ± 1.1 

10 . 2 ± 1.4 

2 . 6 ± 0.2 

c 1.2 

13±3 

23 


c 1.3 


Notes: Data without errors are upper limits; a 2MASS data; b VVV-DR1 data; r No confirmed detection; d VVV-DR1 fluxes in Jy x 10 5 : F ; =761±2 and Fy=3017±4. 
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imising the deviations from the de-reddened flux densities. In 
several cases two models were fitted representing the star at dif¬ 
ferent variability states. Because of the requirement of a very 
red IRAS colour, most stars possess a CSE with a high optical 
depth Tg j^m >10, and only four have 5 < <10. The bolo- 

metric fluxes determined vary by a factor 20 and are in the range 
fbot = 10 -12 -2 x 10 -11 Watt nr 2 . 

The results are given in Table[3] It lists the object name in the 
first column. The best-fit optical depth Tgj llm , and the bolometric 
flux fb 0 i delivered by DUSTY are given in the two following 
columns. The last two columns give luminosities and mass-loss 
rates discussed in the next Section. 

The uncertainties of log fb 0 i are approximately 0.2 dex, 
which we attribute mostly to variability. For the cases where 
models for different variability states are given, the difference is 
~ 0.3 dex. We conclude therefore that the peak-to-peak luminos¬ 
ity variation of the large-amplitude variab le objects is > 2. This 
number is probably a lower limit, as ISuhl (120041) determined for 
two extremely obscured OH/1R stars in the Galactic disk lumi¬ 
nosity variations by a factor 3-3.6. The bolometric fluxes deter¬ 
mined here cannot be individually corrected for variability, but 
in general we expect them to be close to the mean bolometric 
fluxes. The deviations of the observed bolometric fluxes from 
the mean fluxes is at most half of the peak-to-peak variations, 
and very likely even lower in those cases where average bolo¬ 
metric fluxes are available from fits at different variability states. 

The de-reddened SEDs and the models are shown in Fig. 7. 
The observed SEDs are the result of the emission from the 
circumstellar dust with high optical depths, which in addition 
are reddened further by interstellar extinction. The removal of 
the interstellar extinction makes the SEDs bluer and brighter 
than the observed SEDs. As an example, the influence of ex¬ 
tinction on photometry and model fit is shown in Fig. |T| for 
IRAS 17392-3020. We found that the increase in logarithmic 
bolometric flux surpasses 0.2 dex only if A k s > 1.0 mag, which 
applies for four sources in the sample. Therefore, the largest con¬ 
tributions to the uncertainties of the fboi determination is in gen¬ 
eral from the variability and not from interstellar extinction. 

For several sources the DUSTY models predict less flux than 
observed in the near infrared at A < 3 gm and in the far infrared 
at wavelengths 60<A< 160 gm. These sources are marked in 
TableQ]as ‘NIR-Exc’ and ‘FIR-Exc’. The observed excess emis¬ 
sion makes only a minor contribution to the overall luminosity 
and has been neglected for the luminosity determination. 

The near-infrared excess presented by IRAS 17128-3528 
and IRAS 18092-2508 (see Fig. 7) is evident from the unexpect¬ 
edly high brightness in the J, H, and K filters measured by the 
VVV-DR1 survey. We verified that this is not due to missiden- 
tification or contamination of neighbouring objects. These stars 
may have developed deviations of their dust distribution from 
spherical symmetry or their circumstellar shell may have been 
diluted after the end of the strong AGB mass loss, so that the 
light of the central star could escape to form the excess. They 
might currently be starting the post-AGB evolution and may 
later develop double-peaked SEDs as seen in the PPN candi¬ 
dates IRAS 17276-2846 and IRAS 17504-3312 (cf. Sect.IXTl 
Alternatively, they might belong to binary systems, where the 
excess emission comes from the companion. In any case, addi¬ 
tional photometry is required to verify the near-IR excesses for 
these sources. 

The SEDs of IRAS 17292-2727, IRAS 17367-2722, 
IRAS 17367-3633, and IRAS 17583-3346 (Fig. 7) present a far- 
infrared excess. This far-infrared excess emission can not be at¬ 
tributed to photometric errors or variability since consistent flux 


IRAS 17392-3020 



Fig. 1 . Observational photometry and model fits of 
IRAS 17392-3020. For each photometric observation, the 
uncorrected flux density (lower red symbols) and the de- 
reddened flux density for A^, = 1.8 mag (upper blue symbols) 
are given. The best-fitting models are given as dotted (red) curve 
for the observed photometry and as solid (blue) curve for the 
de-reddened photometry. The de-reddened bolometric flux is a 
factor 2.2 enhanced. Symbols are as in Fig. 7. 

densities have been measured at different times by more than one 
survey. Such excess emission has been seen before in a num¬ 
ber of objects, in which the excess can be followed over a large 
wavelength interval extending up to the submillimeter and mil¬ 
limeter ranges. To account for the excess flux an additional dust 
component of larger cool grains is required to model the SEDs 
of these stars (ISanchez Contreras et al.|[2007l) . 

3.3. Luminosity 

We estimated the luminosity from the bolometric fluxes fboi, as¬ 
suming a common distance to all of the sources in the sample 
equivalent to the distance to the Galactic centre (8.0 + 0.5 kpc). 
The luminosity values obtained are tabulated in the fourth col¬ 
umn of Table[3] In case of two fits, representing a brighter and a 
fainter phase, both luminosities are given in the table. 

To quantify the uncertainties of the luminosity determina¬ 
tion, we take the uncertainties of the distances and of the bolo¬ 
metric flux determination into account. Taking the Galactic cen¬ 
tre distance uncertainty and the estimated size of the bulge 
(~2.8 kpc) into account, the distance dependent uncertainty of 
the luminosity is - 35%. This is of the same order as the un¬ 
certainty of the bolometric flux due to variability and extinction 
(see Sect. im Therefore, luminosities of individual objects are 
only accurate to within a factor of ~ 2. 

Figure [2] shows a histogram of the luminosities in our 
Galactic bulge sample. In the case of two model fits, we used the 
average luminosity. The range of luminosities found goes from 
~ 3000 to ~ 30,000 L e , peaking at about 4500 L 0 . 

As mentioned in Section 2, the minimum luminosity re¬ 
quired by AGB stars to be present in the current sample (e.g. 
F 25 > 7 Jy) depends on the optical depth of the shells. Stars at a 
distance of 8.0 kpc and a SED represented by one of the DUSTY 
models would need for > 20 merely a minimum luminos¬ 
ity of 2000-3000 L 0 . For lower optical depths 5 <Tgj m <20 
this minimum luminosity rises to ~ 8500 L 0 . Luminosities of the 
order of a few thousand solar luminosities in Table [3] are there- 
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Table 3. Result of the DUSTY model fit. 


(1) 

IRAS 

name 

(2) 

m 

(3) 

-log( fbolT 

(4) 

L b 

xlO 3 L 0 

(5) 

M 

xl0~ 5 M 0 yr~‘ 

16582-3059 

15.9 

11.70 

4.0 

1.7 

17030-3053 

10.9 

11.75 

3.5 

1.0 

17107-3330" 

6.2 

11.64 

4.5 

1.0 

17128—3528 c 

41.3 

11.30 

9.9 

8.7 

17151—3642 c 

23.4 

11.25 

11.1 

6.0 

17171—2955 c 

10.9 

11.63 

4.7 

1.4 

17207-3632 

34.2 

11.00 

19.8 

12.8 

17251-2821 

23.3-28.2 

11.44-11.80 

7.2-3.1 

1.7-3.4 

17292-2727 

15.9 

11.05 

17.7 

5.7 

17316-3523 

15.9 

11.08 

16.5 

7.6 

17323-2424 

28.2 

11.67 

4.2 

2.7 

17341-3529" 

5.1-6.2 

11.05-11.30 

17.7-9.9 

2.1 -3.4 

17350-2413" 

9.0-10.9 

11.50-11.95 

6.3-2.2 

0.7-1.7 

17351-3429 

13.2 

11.30 

9.9 

2.4 

17361-2358" 

15.9 

11.57 

5.3 

2.1 

17367-2722 

34.2 

11.60 

5.0 

6.8 

17367-3633 

10.9-13.2 

10.80-11.15 

31.4-14.0 

4.7-9.3 

17368-3515 

28.2 

11.56 

5.5 

3.6 

17392-3020 

41.3 

11.01 

19.4 

11.3 

17418-2713 

50.0-50.0 

10.67-11.02 

42.4-18.9 

18.0-40.4 

17428-2438 

15.9 

11.40 

7.9 

3.2 

17495-2534 

15.9-28.2 

10.70-11.00 

39.6-19.8 

11.1-14.0 

17521-2938 

41.3-50.0 

11.60-11.88 

4.9-2.6 

2.3-3.7 

17545-3056 

50.0 

11.77 

3.4 

3.3 

17545-3317 

15.9 

11.50 

6.3 

2.6 

17583-3346 

5.1 

11.40 

7.9 

1.6 

17584-3147 

34.2 

11.45 

7.0 

4.3 

18019-3121" 

28.2 

11.71 

3.9 

2.5 

18040-2726 

34.2-34.2 

11.35-11.65 

8.9-4.4 

2.8-5.7 

18040-2953" 

19.3-23.3 

11.30-11.65 

9.9-4.4 

2.4-4.6 

18091-2437 

15.9-19.3 

11.50-11.85 

6.3-2.8 

1.2-2.3 

18092-2347 

50.0 

11.19 

12.8 

10.6 

18092-2508 

19.3 

11.35 

8.9 

4.0 

18195-2804 

15.9 

11.55 

5.6 

1.9 

18201-2549" 

5.1 

11.25 

11.2 

2.2 


Notes: a f hol in Wattm 2 ; b Luminosity estimated using the distance of 8.0kpc to the Galactic centre. " v exp = 14kms 1 was assumed. Otherwise 
v„ p is derived from either OH or CO radio data. 


fore found only for high optical depths Tgj pm > 10. Stars with 
CSEs of lower optical depths could enter the sample only with 
larger luminosities. 


3.4. Mass-loss rates 

The calculations of mass-loss rates by DUSTY assume that the 
mass loss is radiatively driven. The model rates and terminal out¬ 
flow velocities are given for stars with a luminosity 10 4 L e , fixed 
gas-to-dust ratio r g d = 200 and dust grain bulk density p s = 3 
gcmT 3 . Using the observed luminosities L (Table|3} and the ex¬ 
pansion velocities v exp listed in Table [T] as proxies for terminal 
outflow velocities, the mass-loss rates were adjusted according 
to the scaling relations given in the DUSTY manual (llvezic et al.l 
1 999 ). For stars without OH maser or CO observations v exp = 14 
kms^ 1 was adopted. The mass-loss rates M obtained, are in the 
range ~ 10 5 -3x 10 4 M 0 yr _1 . The rates for individual stars 
are listed in Table [3] Because of uncertainties inherent to the 


DUSTY code of 30% and the uncertainties of the luminosities, 
the error of individual mass-loss rates is also a factor ~ 2. 

The scaling relations adjust the product r g d ■ p s to accommo¬ 
date the observed expansion velocities. With the DUSTY values 
for the gas-to-dust ratio r g d and dust grain bulk density p s , the 
model outflow velocities are generally only about half as large 
as the observed velocities. To increase the model velocities, the 
product needs to be decreased, as \r exp oc L 0 25 • ( r g d ■ p s Y°' 5 - 
As p s = 3 gcm~ 3 is already considered as a c onservatively low 
density for astronomical silicates dDraine & Leell 19841) . the gas- 
to-dust ratio is probably smaller. For the 27 sources with ob¬ 
served expansion velocities (cf. Table Q}. r g d = 44 + 20 is re¬ 
quired to bring the model and observed expansion velocities into 
agreement. This value corresponds to the lower end of the range 
50 < < 1 80 determined for OH/IR stars in the Galactic disk 

bv lJusttanont et all d2006l) . 

The mass-loss rates in general are rather high. To check 
their reliability, we compared them with those of iGroenewegenl 
(2006;), who calculated mid- and far-infrared colours for mass 
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'Luminosity (xlCTL,;,) 

Fig. 2. Luminosity distribution of the Galactic bulge sample. 
Luminosities were obtained from the model fitting of the SEDs. 
In the case of two fits for the same source, the average has been 
taken. 
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Fig. 3. Mass-loss rates and luminosities derived from the 
DUSTY models fitting of the AGB star SEDs in our Galactic 
bulge sample. In cases where ranges of luminosities and mass- 
los rates are given in Table [3 we plot the average. Stars with 
envelopes of unknown chemistry are plotted with asterisks. The 
classical limit for mass-loss rates caused by radiation pressure 
on dust is shown by a dashed line, and the empirical mass-loss 
rate upper boundary for AGBs_and red supergiants in the Large 
Magellanic Cloud found by Ivan Loon et al.l (119991) is plotted as 
a dotted line. Also the maximum luminosity achievable by AGB 
stars is marked by the vertical dash-dot line. 
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losing AGB and post-AGB st ars ba sed on own models. From 
the different dust compositions iGroenewegeni ( 2006 ) employed, 
we used his pure silicate models for comparison. The range of 
mid-IR colours he considered covers only about 40% (the blue 
part) of the colours of the GB sample, but among these the mass- 
loss rates are in reasonable agreement. The mean ratio between 
Groenewegen’s and our mass-loss rates is 0.9 ±0.4. 

Figure [3] plots the mass-loss rates against the luminosities 
of the sample. Because of scaling relations the two quantities 
are not independent from each other, and in the cla ssica l case 
a relation M r i„ SS ir — L ■ ( cv exp )~ l ccL 0 75 (Ivan Loon et al.l[l999t 
lElitzur & Ivezidl200ll) is expected. The scatter of the mass-loss 


rates for a given luminosity is caused mainly by the scatter of 
optical depths Tgj required to model the SEDs and to a lesser 
degree by the variations of the expansion velocities. In general, 
the luminosities and mass-loss rates cover the range between the 
classical limit in the (M, L)-plane and the empirical limit found 
bv lvan Loon et al.l(ll999l) in the Large Magellanic Cloud (LMC). 
The classical limit is determined by the assumption that the pho¬ 
tons transfer momentum to the dust by a single scattering event. 
This limit is superseded however at high optical depths, where 
photons experience multiple scattering events, before escaping 
from the dust shell. Compliance with this empirical limit was 
found also for the AGB stars in the central regions of M33, 
which have_a near solar metallicity comparable with the Galactic 
bulge (iJavadi et al .1120131) . 

Few stars exceed the empirical limit. For example, 
IRAS 17418-2713 has the largest mass-loss rate in the sample 
of 3 x 10 4 M G yr , and exceeds the empirical limit by a factor 
~ 1.5. The very steep rise of the SED in the near-infrared re¬ 
quires a high optical depth for the model dust shell, which leads 
to a high mass-loss rate in combination with the inferred high 
luminosity of ~ 30,000 L G . However, given the uncertainties for 
the luminosities and mass-loss rates, a systematic difference of 
the upper limit of the (M, L)-relation of AGB stars in the bulge 
and the LMC cannot be inferred. 


4. The nature of the sample 

4.1. The distribution of luminosities and mass-loss rates 

Our sample is made up of the most reddened AGB stars of the 
Galactic bulge known so far. We found a large range of lumi¬ 
nosities with a peak in the luminosity distribution at L ~ 4500 L G 
consistent with low-mass progenitors, and with a tail extending 
to high luminosities. 

This result is in good agreement with other studies of less 
reddened AGB stars in the Galacti c bulge. The luminosity dis¬ 
tribution of lBlommaert et all ( 1998 ) for a sample of OH/IR stars 
located close to the Galactic centre spans a similar range and 
has a median value of ~4600 L G . The median l umino sity of 
Galactic centre OH/IR stars monitored by I Wood et ail (119981 ) 
for variability is ~ 6400 L G , and their range of luminosities ex¬ 
tends to L > 30,000 L g consistent with an upper limit for the 
mass range > 6 M G . Both studies relied on measurements in the 
near-infrared (K, L-band), making the colours of their stars bluer 
than in our sample. 

Intermediate in colours between the Galactic centre OH/IR 
stars of the just above mentioned studies and the colours of 
our samp le are the OH/IR stars observed by the ISOGAL sur¬ 
vey dOrtiz et al.1120021) and of the MSX-selected selected sam¬ 
ple of AGB stars (lOiha et al .1120071) . Their luminosity distribu¬ 
tions peak around 8000 L G , and objects with higher luminosi¬ 
ties are found in both sam ples. The red IRAS sources studied 
by Ivan der Veen & Habind dl990h are very similar to our sam¬ 
ple. Their sources cover the colour range 0.0 < [12] - [25] < 1.3, 
which overlaps in its red part with the colour range ([12] - 
[25] > 0.75 mag) of the present sample. Most of their sources are 
bluer than ours, with only seven sources in common. They find 
the luminosities to be strongly peaking at 5000-5500 L 0 and ob¬ 
served a tail of high luminosity sources extending to well above 
20,000 L 0 . 

We therefore conclude that the luminosity distributions of 
the different samples of AGB stars observed in the Galactic 
bulge are consistent with our result in showing that the majority 
of the stars have luminosities well below 10,000 L 0 , but that stars 
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with luminosities up to the AGB limit are present as well. There 
is no evidence for differences in the luminosity distribution with 
colour, which would indicate the presence of mass segregation. 

Mass-loss rates in the range from 10 7 to 1 0 4 Mp vr 1 for 
Galactic bulge A GB stars were derived by rOiha et all (120071) . us¬ 
ing the models of lGroenewegenl d2006l) . There are nine sources 
in common with our sample. The deviations in the mass-loss 
rates are a factor 0.6-1.5, except for IRAS 17418-2713, in 
which the rates differ by a factor of 10. lOiha et aTI (l2007h un¬ 
derestimate the rate in this case, because of their adoption of 
a mean luminosity for all objects to scale the model mass-loss 
rates. IRAS 17418-2713 is actually a factor ~4 more luminous, 
which explains at least part of the difference in the mass-loss 
rates. 


4.2. Predictions from AGB evolutionary models 


High mass-loss rates of the order 10~ 5 Moyr” 1 or higher 
(dubbed superwind hereafter) imply that the stars lose most of 
their mass in a relatively short time. Actually the mass-loss rate 
determines the stellar lifetime on the AGB. According to the 
models, during the evolution on the TP-AGB, lasting ~ 0.25 - 
2.2 million years, luminosities and mass-loss rates steadily in¬ 
crease, until the last ~ 60 ,000-120,000 years when the stars en¬ 
ter the superwind phase (tVassiliadis & Woodlll993l) . These time 
ranges vary slightly depending on the metallicity assumed and 
how the mass-loss process is implemented in model calcula- 

tions dBlockerll 19951: IMarigo & Girardil2007tlWeiss & Fergusonl 

l2009h . The stars of our sample have mass-loss rates in the range 
from (~l-30)xl0~ 5 M 0 yr 1 , characteristic of the superwind 
phase, and we assume that they will leave the AGB within the 
next ~ 10 5 years. Hence, a comparison with stellar evolutionary 
models has to assign these stars to the low-temperature, high- 
luminosity ends of the evolutionary tracks on the AGB for low- 
and intermediate mass stars. 


The final phases of these tracks on the AGB are characterised 
by similar effective temperatures of 2500-3000 K, and mod¬ 
els for different zero-age main-sequence (ZAMS) masses are 
then basically distinguished by the fina l luminosity reached and 
in part by their chemica l comp osition (IMarigo & Gira rdi 2007; 
Weiss & Fergusonl2009l) . While lVassiliadis & Woodl l 19931) and 


Blockei ( 1995b do not discuss the conversion of O-rich to C- 


rich chemistry in their A GB models due to the third dredge-up, 
IMarigo & Girardil (120071) predict that stars in the mass range 
~2. 0-4.0 M g with solar metallicity (Z = 0.019), or ~1.5- 
4.0 M 0 with typical metallicity for th e LMC, (Z = 0.008), end 
as C stars. In I Weiss & Fergusonl (120091) , the corresponding mass 
range for solar metallicity is 1.8-3.0 M 0 and, for Z = 0.008, no 
model with Mms <5.0 M 0 ends O-rich. For most stars of our 
sample (80%), however, the O-rich chemistry is evident from the 
presence of masers and/or a signature of silicates in the infrared. 
Hence, if these stars are indeed close to the end of their TP-AGB 
evolution, they should have evolved either from intermediate- 
mass stars with Mms >4.0 M Q experiencing HBB, or from low- 
mass stars of at least solar metallicity. This conclusion would be 
avoided by the suppression of carbon star formation for the rest 
of the masses, which may occur for higher than solar metallicity 
or for increased helium content (iKarakasll2014l) . 


Figure [4] compares the locations of the GB sample 
stars in the Hertzsprung- Russell diagram with the models of 
iBertelli et al.l (120081. [2009]) , describing the TP-AGB stellar evo¬ 
lution for ZAMS masses 0.7 < Mms <6.0 M e , assuming so¬ 
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Fig. 4. Location of the GB sample stars in the Hertzsprung- 
Russ ell diagram overlayed on e volutionary tracks on the AGB 
from IBertelli et al.l (120081120091) for solar composition. Tracks 
are for main-sequence masses 0.7 < Mms ^ 6.0 M G . Segments of 
the tracks, where the photospheric chemical composition is C- 
rich, are shown in red (dotted). Stars with envelopes of unknown 
chemistry are plotted with asterisks. Average luminosities of the 
GB sample stars are taken from Tabled For the unknown effec¬ 
tive temperatures random values between 2500 and 2800 K were 
adopted. 


lar composition (Z = 0.017, Y = 0.260 The effective temper¬ 
atures T e ff of the stars are not known, but we assume that 
they are < 3000 K. To avoid crowding in the figure, we as¬ 
signed random effective temperatures to the stars in the range 
2500 < T e ff <2800 K. Contrary to the model predictions, the 
mass range 1.1^ Mms ~ 6.0 M 0 is continuously covered by the 
sample, including the range in which the stars should end as C- 
rich. Most of the stars, in which the chemistry of the envelope 
is unknown, have relatively low luminosity in line with our as¬ 
sumption that they are all oxygen rich. 

4.3. The high-luminosity group 

The luminosity distribution of the GB sample (Fig. 0 has a 
steep increase at low luminosities, peaks at ~4500 L 0 , and 
falls off less steeply towards higher luminosities. For the fol¬ 
lowing discussion, we name sources with L> 10,000 L 0 the 
‘high-luminosity group’ and those with L < 7000 L 0 the ‘low- 
luminosity group’. 

In the ‘high-luminosity group’ two sources 
(IRAS 17418-2713 and IRAS 17495-2534) have the highest 
luminosities (~ 30,000+ 10,000 L 0 ; Table |3j. Apart from the 
possibility that they could be foreground sources and hence 
are actually less luminous, their luminosities are compatible 
with intermediate-mass stars ( Mms ~ 4 - 6 M 0 ) experiencing 
HBB, preventing them from becoming C stars. Both show the 
absorption feature from amorphous silicates at 10/rm (Table 
|T}. confirming their O-rich chemical composition. Furthermore, 
IRAS 17418-2713 sho ws cry stalline silicate emissi on fea¬ 
tures in the infrared (iGarcla-Hernandez et al.l 12007 ). while 
IRAS 17495-2534 is the only sourc e known so far showing 
these features in absorption dSneck et al.l 120081) . Speck et al. 
argue that the rarity of such a distinctive feature is due to their 

5 Available at http://stev.oapd.inaf.it/YZVAR/ 
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emergence only in massive AGB stars, which is consistent with 
our findings. 

The rest of the group, making up about one quarter of the 
whole GB sample, has luminosities between 10,000-25,000 L 0 . 
These luminosities are predicted for the mass range, where for 
solar composition, stars end as C-rich. For lower metallicities, 
this group will even increase, as conversion to C-rich chem¬ 
istry also occurs for smaller masses. However, masers and sil¬ 
icate features in the infrared are observed in most stars from 
this group, ruling out a C-rich chemistry. It is unlikely that a C- 
rich central star coexists with OH maser emission in all these 
stars, because such a configuration can last only for ~ 1000 
years, which is the time needed to replace the O-rich material 
in the OH maser shell at Roh ~ 5 xlO 16 cm, assuming typical 
outfow velocities of 15 kms -1 . As the effective temperatures 
are unknown, these stars in principle could have temperatures 
T e ff > 3000 K, assigning them to an earlier evolutionary phase, 
where the chemical composition is still O-rich, and the lumi¬ 
nosities are almost t he same (cf. F ig. HJi. H owever, none of the 
model calculations (Vassil iadis & Woodl 19 93t iBlockeil 1 1995b 


iMarigo & Girardil2007 : Weiss & Fergusonl2009t) predict the ob¬ 

served high mass-loss rates at less advanced phases of TP-AGB 
evolution. 

The presence of O-rich objects at the tip of the AGB with 
luminosities, corresponding to masses where the stars should be 
already converted to C stars, may have two possibility. Either 
phases of very high mass-rates can occur, before the O —»C 
transition is experienced, or for some of the stars from this mass 
range the transition might be delayed or avoided. The first pos¬ 
sibility is unlikely except in the case of higher mass stars, be¬ 
cause it is likely that very little mass is left over after a su¬ 
perwind phase to continue the evolution with further thermal 
pulses, allowing an O —> C transition before the stars leave the 
AGB. Th erefore, the second possib ility may be more rewarding 
to follow. I Weiss & Ferguson ( 1200 91 find that models with super¬ 
solar metallicities can avoid the conversion to C/O > 1 before 
the end of the AGB simply because they start with smaller C/O 
ratios, and therefor e m ore thermal pulses are needed to reach 
C/O ~ 1. Also, iKarakasi (120141) finds that at a super-solar metal- 
licity Z=0.03, no carbon stars are formed except in a narrow 
mass range of 3.25-4 M 0 . Therefore, in environments where a 
mix of metallicities is present, a mix of stars originating from 
similar ZAMS masses, but with different chemistries, might also 
be present. It would imply that most of the objects from the GB 
sample were born in bulge regions with higher than sola r me tal¬ 
l icities . A similar line of argument was forwarded bv lJura et alj 
( 1199 31) to explain the differences in the period distribution of 
Mira variables in the solar neighbourhood and the Galactic cen¬ 
tre. They suggested that Mira variables have higher metallicity 
in the Galactic centre and are therefore oxygen rich, while stars 
of the same mass in the solar neighbourhood are less metal rich 
and already converted into carbon stars. 

The metallicity argument can be tested for consistency in 
comparison with the spatial distribution of the GB sample with 
the global metallicity distribution in the Galact ic bu lge. Fig. 
[5] shows the metallicity map of Gonzalez et al.l ( 2013 1). which 
clearly confirms the vertical metallicity gradient with metal-rich 
stars dominating the inner bulge in regions closer to the Galactic 
plane (| b \< 5°). The GB sample members are overplotted. They 
are completely confined to the more metal rich part of the bulge. 
Although the mean metallicities in the 30'x45' bins of the metal¬ 
licity map do not reach super-solar metallicities, higher metal¬ 
licities may occur on smaller spatial scales, where the GB sam¬ 
ple members may be located. Therefore, at least this comparison 
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Fig. 5. Distribution of the GB sample stars overlayed on the 
d istributi on of mean metallicities in the Galactic bulge of 
iGonzalez et aL (120131) . The different symbol sizes mark the 
members of the low- and high-luminosity group and the stars 
with luminosities in between. 


does not contradict the model predictions that the chemistry of 
intermediate-mass stars at the termination of TP-AGB evolution 
depends on the metal content of the environment in which they 
were born. 

4.4. The low-luminosity group 

Almost half of the sample has luminosities L<7000 L G 
(logL<3.85) and mass-loss rates (l-5)xl0“ 5 M 0 yr -1 . These 
kinds of combinations of relatively low luminosities and very 
high mass-loss rates are not predicted by the models either, at 
least as long the stars are considered to be on the TP-AGB and 
in between two thermal pulses (the interpulse phase). Assuming 
that the stars are in the final phase of their AGB evolution (the 
last couple of interpulse periods) and making the requirement 
that the luminosities do not surpass 7000 L G , these stars must 
have descended from 1.1 ^Mms ~ 1 -8 M 0 main-sequence stars 
(Fig.[4]i. 

In models of TP-AGB evolution, mass-loss rates exceeding 
10- 7 M g are not assumed to occur before large amplitude pulsa¬ 
tions of the stellar envelope have been established. Most mem¬ 
bers of the GB sample have not been monitored for variab ility, 
except the few sources in common with Ivan der Veen & Habind 
( 1990 ). for which they report large amplitude variability with 
periods P> 750 days. The other sources usually have high IRAS 
variability index (> 95, see Table [TJ- indicating that they are 
also, in the majority, large amplitude variable stars. We there¬ 
fore exclude the remote possibility that the members of the low- 
luminosity group are intermediate-mass stars ( Mms i 3-0 M G ) on 
the early AGB preceding the TP-AGB, although the observed 
lumino sities (L<7 000 Lo) are achieved by models during this 
phase dBertelli et al.ll2008» . 

During the TP-AGB, the mass-loss rates steadily increase, 
modulated by the intervening thermal pulses, which affect the 
surface luminosity and radius of the stars. During the inter¬ 
pulse phase, the stellar luminosity reaches its maximum shortly 
before the next thermal pulse happens and is_ referred here 
to as q uiescent luminosity. The models of IVassi liadis & Woodl 
( 1993 ) provide the mass-loss rates (> 10 5 M G yr 1 ) observed 
for Mms > 2.0 M 0 during the last couple of interpulse peri¬ 
ods, but the correspondent model luminosities are too high. 
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[ 3 . 6 ] [ 8 . 0 ] 


Fig. 6. Location of the GB sample stars (red filled circles) in the 

[3.6] —[8.0] vs. [8.0] colour-magnitude diagram assuming a dis¬ 
tance corresponding to the LMC. The boundaries of the locations 
of the different type of AGB stars (O-rich, C-rich and extre me) a s 
defined by the SAGE project are taken from lBlum et al . (12 0061) . 
In addition, the ‘evolved s tars’ in the LMC (plus signs) identified 
by iGruendl & Chul (120091). and the spectrosc opically confirmed 
C stars (asterisks) from lGruendl et all (120081) are overplotted. 

Blocked d 19951) presents models with Mms = E 3, 4, 5, and 
7 M q . While the 1 M G model does not reach the TP-AGB 
phase, the Mms = 3 M e model has quiescent luminosities in the 
required range, but does not reach mass-loss rates surpassing 
10~ 5 M 0 yr _1 . Models for larger masses do reach the required 
mass-l oss rates, but are again too lu minous. iMarigo & Girardil 
(120071) and lWeiss & Fergusonl (120091) also arrive at the same re¬ 
sult: the high observed mass-loss rates are not predicted for the 
mass ranges corresponding to the observed quiescent luminosi¬ 
ties. Given the short times the therma l pulses last (~ several hun¬ 
dred vears. fVassiliadis & Woodll993l) . it is unlikely that all these 
stars are currently experiencing a boost of the mass-loss rate fol¬ 
lowing the brief increases of the luminosity. 

A possible explanation for the discrepancies between the ob¬ 
served properties and the model predictions may originate in the 
assumptions concerning the relation between mass-loss rates and 
other properties in the TP-AGB models. As no complete mass- 
loss theory for AGB stars is available, mass-loss rates are empiri¬ 
cally coupled to pulsational period. The periods are theoretically 
determined from relations involving mass and radius, which are 
provided by the models. While these relations provide resonable 
values in general, there might be brief phases where the pre¬ 
dictions are insufficient. At least for the low luminosity part of 
our sample, which had main-sequence masses Mms < 2 M G , the 
models predict mass-loss rates that are too low in the brief final 
evolutionary phase. 

4.5. Comparison with the AGB population in the LMC 

In the LMC, the AGB population has_been studied extensively 
by the SAGE project dBlum et al.l 12006 ). The chemistry of the 
population was determined from photometry using basically the 
J— [3.6] vs. [3.6] colour-magnitude diagram (CMD). In Fig. [6] 
the boundaries for different chemistries are delineated in the 

[3.6] -[8.0] vs. [8.0] CMD, which is more appropriate for very 
obscured stars. While C stars are in general found to be red¬ 
der than O-rich stars, Blum et al. also define a group of ‘ex¬ 


treme stars’ with even redder colours (J - [3.6] >3.1 mag), for 
which the chemistry could not be determined. The few known 
stars of this group, which have been observed spectroscopically, 
are predominantly C-rich. Along with the model predictions that 
the lowest effective temperatures and hence the largest mass- 
loss rates are achieved when the stars become C-rich, the ‘ex- 
treme stars’ are uniformly classified as p redominantly C-rich 
(ISrinivasan et al.l 1201 It IRiebel et alJl2012l) . However, a minor¬ 
ity population of O-rich stars is not excluded and contamination 
with re d supergiants and young stellar objects i s probably also 
present (ISrinivasan et al.ll2009HBlum et al.ll2014l) . 

We determined synthetic apparant magnitudes in the 
3.6] and [8.0]_bands of the GLIMPSE photometric system 
Churchwell et al.l l2009t) for all stars of the GB sample using 
the DUSTY model fits. Comparison with the actually observed 
magnitudes by GLIMPSE gives an uncertainty of ~ 1 mag. To 
comp are the location of the sample in the CMD of iBlum et al.l 
(120061) . a distance modulus of 4.0 mag was used to account for 
the distance difference between the LMC and the Galactic bulge. 
As shown in Fig.[6j the GB sample stars are located in a sparsely 
populated region extending the branch populated by the ‘ex¬ 
treme stars’ towards redder colours [3.6 1 - [ 8.01 > 2 .5 mag. This 
region of the CMD was also found by IGruendl & Chul d2009l) to 
be populated by ‘evolved stars’. The reddest of them (dubbed 
EROS for ‘Extremely Red Objects’) were spectroscopically ob¬ 
served by IGruendl et al] (120081) . showing C-rich chemistry. It 
follows that while the SAGE ‘extreme sta rs’ and the even redder 
‘evolved stars’ of IGruendl & Chul d2009l) in the LMC are pre¬ 
dominantly C-rich, the corresponding population in the Galactic 
bulge represented by the GB sample seem to be O-rich. 

IRiebel et al .. d2012t) determined dust mass-loss rates for the 
SAGE stars using model SEDs. Almost all ‘extreme stars’ were 
matched with C-rich models. Most dust mass-loss rates for O- 
rich stars are in the range 10~ n - 10 9 M G yr 1 , and for C-rich 
stars > 10~ 10 M 0 yr _1 . With a gas-to-dust ratio of ~200 their 
total mass-loss rates are well below the rates achieved during 
the superwind. The separation of O-rich and C-rich sources in 
the SAGE CMDs is very likely valid only as long as mass-loss 
rates are low, the CSEs are optically thin and therefore the stel¬ 
lar temperature is able to influence the infrared colours. As soon 
as the mass-loss rates increase, and especially during superwind 
phases, the stars migrate in CMDs towards redder colours, inde¬ 
pendent of their current chemistry. IGruendl et al. ( 2008 !) find a 
mean luminosity of 7100 L 0 for their EROS sources, and con¬ 
clude that these stars must have descended from 1.5-2.5 M G 
main-sequence stars. The total mass-loss rates derived are in the 
range of (4-23) x 10 5 M e ,yr 1 , much larger than derived for 
the SAGE AGB stars. The EROS luminosities, main-sequence 
masses, and mass-loss rates are in perfect agreement with those 
derived here for the GB sample. The main difference is the chem¬ 
ical c ompositio n. As_we previously concluded for the GB sam¬ 
ple, IGruendl et al.l (2008) noted that the mass-loss rates of the 
EROS sources typically exceed the maximum expected for both 
O-rich and C-rich AGB stars with luminosities in the observed 
range. 

5. Conclusions 

The reddest bright 25 pm sources towards the Galactic bulge 
with IRAS colours expected for AGB stars were found to span 
a luminosity range from ~ 3000 to ~ 30,000 L e and mass-loss 
rates between ~ 10 5 - 3 x 10 4 M 0 yr 1 . For the lower lumi¬ 
nosity group with L<7000 L 0 the main-sequence progenitor 
masses are in the range 1.1-1.8 M 0 , assuming solar metallic- 
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ity. With mass-loss rates (l-5)x 10 5 M 0 yr their stellar en¬ 
velopes are lost in mere 10,000-200,000 years and AGB evo¬ 
lution is terminated. One concludes that the heavily obscured 
AGB stars do not necessarily descend from massive AGB stars 
(Mms >4M q ) only. At higher luminosities (10,000-25,000 L 0 ), 
AGB evolutionary models predict that for metallicities similar to 
the Sun or less these high-mass loss stars should already be con¬ 
verted to C-rich chemistry. On the contrary, in the bulge these 
stars are mostly O-rich. The models predict that the conversion 
to C-rich chemistry can be delayed or inhibited for higher than 
solar metallicities. Our sources might have been born in regions 
of the bulge, where several generations of stars have leftover pro¬ 
cessed matter with higher than solar metallicity. Their spatial 
distribution favourably coincides with regions of higher mean 
metallicities in the Galactic bulge. 
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IRAS 16582-3059 


IRAS 17128-3528 




IRAS 17030-3053 


IRAS 17151-3642 




IRAS 17107-3330 


IRAS 17171-2955 




Fig. 7. Observation and model fits after correcting for interstel¬ 
lar extinction. 2MASS: circles; VVV-DR1: stars; WISE: trian¬ 
gles pointing up; GLIMPSE: sails; MSX: squares; AKARI: di¬ 
amonds; IRAS: crosses; Upper limits: filled triangles pointing 
down; Unconfirmed detection: open triangles pointing down; 
Model: solid line (in case of two models, the solid line represents 
the brighter and the pointed line the fainter phase (see Section 

ED). 
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Fig. 7. continued. 
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Fig. 7. continued. 
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Fig. 7. continued. 
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Fig. 7. continued. 


Fig. 7. continued. 
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Fig. 7. continued. 
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Fig. 7. continued. 
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